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ABSTRACT 

On April 28-29 2010, the Solar Terrestrial Relations Observatory (STEREO)- 
B/Extreme Ultraviolet Imager (EUVI) observed four homologous large-scale 
coronal waves, so-called 'EIT- waves', within eight hours. All waves emerged from 
the same source active region, were accompanied by weak flares and faint CMEs, 
and propagated into the same direction at constant velocities in the range of 
~ 220 — 340 km s~^. The last of these four coronal wave events was the strongest 
and fastest, with a velocity of 337±31 kms~^ and a peak perturbation amplitude 
of ~ 1.24, corresponding to a magnetosonic Mach number of Mms ~ 1.09. The 
magnetosonic Mach numbers and velocities of the four waves are distinctly cor- 
related, suggestive of the nonlinear fast-mode magnetosonic wave nature of the 
events. We also found a correlation between the magnetic energy build-up times 
and the velocity and magnetosonic Mach number. 



-2- 



Subject headings: Sun: corona — Sun: coronal mass ejections (CMEs) — Sun: 
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Introduction 



Large-scale propagating disturbances in the solar corona have been subject of extensive 
studies for more than twelve years. For the first time, these events we re observed by the 
Extreme-ultraviolet Imaging Telescope (EIT; Delaboudiniere et al. Il995) onboard the Solar 
and Heliospheric Observatory (SoHO; iMoses et al.l 119971 : [Thompson et al.lll998[ ). and are 
thus commonly called 'EIT- waves'. They appear as diffuse coronal brightenings, forming 
circular wave fronts, traveling hundreds of Mm without hindrance under quiet-Sun condi- 
tions. 

At present there are two competing groups of models based on completely different 
ideas of the physical nature of EIT, or more generally, EUV waves. The first theory de- 
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Veronig et al.l l2006l ). Several characteristics of coronal waves can be explained by 



this theory: (1) the propagation perpendicular to magne tic field lines (e.g. 
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Warmuth 
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1998h. (2) the pulse broadening and amplitude drop off (jWills-Daveyl l2006l : 
Veronig et al.l l2010r). ( 3) the reflection and refractiori at regions of high Alfven ve lo city 
jThompson et al.lll999l : fveronig et al1l2008l : Eong et al.ll2008[ lOopalswamv et al.ll2009h . The 
second group of models, however, consider these coronal bright fronts as no waves at all, 
but as the signature of a large-scale rnagnetic reconfiguration of field lin es during a CME 
lift-off (e.g. belannee fc AulanierlE999l : Ichen et al.ll2002[ lAttrill et allboO?! ) . In these pseudo- 
waves models, coronal waves could exclusively occur in conjunction with CMEs. This as- 
sumption is supported by statist ical studies, showing a close relation between waves and 
CMEs (e.g. Biesecker et aP 2002[). These mod els can also explain stationary brightenings 



( jPelannee fc Aulaniei 
coronal bright fronts 



19991: 
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20051 ) . For recent rev iews we refer to 



The limitations of EIT observations, especially the low imaging cadence of ~12-15 
minutes, were a major contributor to the difficulties d etermining the phys ical nature of 
coronal waves. The Extreme Ultraviolet Imager (EUVI; IWuelser et al.l 120041 ) instruments, 
which are part of the Sun Earth Connection Coronal and Heliospheric Investigation (SEC- 
CHI; [Howard et al.l |2008|) suite onboard the twin Solar-TErrestrial Relations Observatory 
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(STEREO; iKaiser et al.l 120081 ) spacecraft overcome these limitations. EUVI observes the 
entire solar disk as well as the corona up to 1.7 Rq in four spectral channels (He II 304 A: 
T ~ 0.07 MK; Fe IX 171 A: T ~ 1 MK; Fe XII 195 A: T ~ 1.5 MK; Fe XV 284 A: T ~ 2.25 
MK). The high imaging cadence down to ~ 75 seconds, the large FoV as well as the ability 
of simultaneous observations from two vantage points, provide us w ith new insights into 



the three-dimensional kinematics and dynarnics o f EUV waves (e.g. iKienreich et al.l 12009 
Patsourakos fc Vourhdasll2009l : IVeronig et al.ll2010l ). 



In this letter we present the first case of homologous EIT waves ever reported in the 
nearly 14 years of continuous studies of large-scale coronal waves. The four homologous 
waves, observed by STEREO-B within a period of eight hours, were launched from the same 
AR, propagated in the same direction, and their wave fronts were similar in both shape and 
angular extent. These events allow us to perform a detailed study of the physical charac- 
teristics of large-scale coronal waves, as these homologous waves are initiated and propagate 
under similar coronal background conditions. Hence, our analysis of the correlation between 
the wave pulse parameters is not influenced by the limiting factor of differing coronal plasma 
conditions, persistent in previous comparative studies of EIT waves. 



2. Data 



The events under study occurred on 2010 April 28-29, and were observed by STEREO- 
B (henceforth ST-B), positioned 70.5° behind Earth on its orbit around the Sun. For our 
analysis, only the high-cadence EUVI imagery in 195 A (cadence of 5 minutes) was suitable, 
since the 171 A and 284 A filtergrams were recorded only every 2 hours, and the 304 A data 
revealed no wave signatures at all. In additi on, the white-light co ronagraph (WLC) images 
of the CORl-B and C0R2-B instruments (iHoward et al.l 120081) and of the Large Angle 
and Spectrometric COronagraph (LASCO; iBrueckner et al.l Il995l ) C2 instrument onboard 
SOHO were used to identify signatures of associated CMEs. The EUVI 195 A filtergrams 
were reduced using the SECCHLPREP routines available within Solarsoft. Additionally 
we differentially rotated the images to the same reference time. In order to enhance faint 
coronal wave signatures, we derived running ratio (RR) images dividing each direct image 
by a frame taken 10 minutes earlier. 
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3. Results 

3.1. Wave characteristics 

The four EUV waves under study were launched on 2010 April 28-29 within a period 
of 8 hours (see Fig. [1]). Movie 1 shows ST-B 195 A direct and running ratio full disk images 
covering the full observation period. Figure 1 consists of snapshots from the movie illustrating 
the evolution of all four waves. In each column we present two fronts of each event, with the 
upper panels showing the wave front with the maximum intensity amplitude, and the lower 
panels showing the propagating wave front 20 minutes later. All four coronal disturbances 
are launched from the same AR NOAA 11067, propagate in the same direction and have a 
similar appearance and angular extent, thus it is appropriate to call them homologous. The 
time interval between the onsets of the successive waves steadily increased: ~1.75 hr (wave 
1-2), ~2.5 hr (wave 2-3) and ~3.5 hr (wave 3-4). In order to obtain better insight into 
the onset of the coronal waves, we studied a small area [FoV: 600" x 600"] centered on the 
source AR. Figure [2]^f,g) gives one snapshot of this region during the launch of wave 4, the 
strongest one of the four events. During the full observation period we recorded persisting 
dynamical processes in the extended loop system of NOAA 11067 (see movie 2 available as 
online supplement). 

The resemblance of the events suggests that the background coronal field has not 
changed noticeably within these eight hours. Moreover, this implies that the corona is 
disturbed by the EIT waves only for a short period of time and returns to its pre-event 
condition on a timescale of tens of minutes. In at least two cases, waves 3 and 4, there is 
evidence that the main perturbation was followed by an additional weak disturbance about 
20 minutes later (see movie 1). 



3.2. Wave kinematics and perturbation characteristics 

We analyzed the kinematics of the waves, each treated as unique event, employing two 
different methods. Firstly, we visually determined the wave fronts in a series of RR images. 
Secondly, we examined the perturbation amplitudes of the wave. In both methods we focused 
on a 60° sector on the solar surface, where all four waves were distinctly observable (see Fig. [1] 
panel 7; yellow curves). 

In the first method the wave fronts were tracked manually, and their center was ob- 
tained from a circu lar fit to the earliest wave fronts carried out in spherical coordinates (see 



Veronig et al.ll2006l ). The mean distance of each wave front from the thus determined center 
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is calculated along the solar surface. In the second method we subdivided the solar surface 
into spherical segments of equal width concentric around the wave center obtained by method 
one. Plotting the average intensities versus mean distances of all segments gives one intensity 



profile per RR image (cf. iMuhr et al.ll2010[ ). In these perturbation profiles, shown in Fig. H] 
(top panels), the wave front presents itself as distinct bump above the background intensity. 
In each case, the perturbation amplitude reaches its maximum ~5 minutes after the onset 
of the wave (see Fig. H]). As example, the evolution of the perturbation amplitude of wave 
4 is plotted in Fig. [3] (bottom panel). From the wave perturbation profiles we extracted the 
foremost position of the wave front, defined as the point at which the gaussian fit to the 
profile falls below the value of ///q = 1.02 (blue dotted line in top panel of Fig. H]). We note 
that the visually tracked distances match well the position of the wave front obtained from 
the perturbation profiles (Fig. H]). 

In Fig. |3] (top panel) we show the kinematics for wave 4 derived with both methods 
together with error bars. Linear as well as quadratic least squares fits were applied to 
the time-distance data, both yielding similar velocities of ~ 340kms~^. The parabolic 
fit suggests a small deceleration of — 10ms~^. In order to distinguish whether the linear 
or parabolic fit better represents the data, we derived the 95% confidence interval for the 
linear fit. The quadratic fit lies within the error bars and the confidence interval, thus it is 
reasonable to represent the wave's kinematics by the linear fit with a constant velocity of 
337 ± 31 kms^^ over the full propagation distance up to 800 Mm. 

We found similar results for the other three wave events under study, which is illustrated 
in Fig. m (bottom panels) . The velocities Vc of the four homologous waves are: 257±22 km s^^ 
(wave 1), 219 ±18 km s"^ (wave 2), 249 ±18 km s"^ (wave 3) and 337 ±31 km s"^ (wave 4). 
All peak perturbation profiles together with a gaussian fit are shown in Fig. H] (top panel). 
The perturbation profiles of all four waves are steepening and show an increase in amplitude 
in the early phase of their evolution until the peak perturbation amplitude Amax is reached. 
The values Amax of all four waves are: 1.15 (wave 1), 1.1 (wave 2), 1.14 (wave 3) and 1.24 
(wave 4). 

The bright fronts of coronal EUV waves are in general caused by a local temperature and 
density enhancement (plasma compression). Assuming that the change in density is more 
important than that in temperature, one can derive an estimate of the density jump from 
the intensity amplitude A = I/Iq, N/Nq ~ (///q)^/^. This implies for the peak amplitudes 
of the four homologous waves a maximum density jump of = N/Nq = 1.07 (wave 1), 1.05 
(wave 2), 1.07 (wave 3), and 1.11 (wave 4). 
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3.3. Associated CMEs and flares 

The analysis of EUVI, CORl and C0R2 images, covering both days of April 28-29, 
revealed no clear evidence of associated CMEs. However, LASCO C2 recorded four faint 
CMEs, each entering the C2 FoV about 45 minutes after the first observation of a wave front 
by EUVI (cf. Fig. [2] (a)-(d) and movie 3). In the LASCO CME catalog just CMEs 2-4 
are listed and classified as poor, only visible in C2 at position angles ~85°. CMEs 3 and 
4 reveal additional bright features following the actual CME leading edge with at least the 
same speed. These features resemble rather small scale ejecta than prominences trailing the 
leading edge and could be related to the above mentioned weak disturbances following wave 
3 and 4 (see movie 1). As the center of the source AR lies ~15° behind the limb, as seen 
from SoHO, we experience low projection effects in the derived CME plane-of-sky velocities, 
which lie in the range of 140 — 190 km s"^ (see Fig. |2] (e)). Setting the onset site to the solar 
limb and onset time to the moment, EUVI observed the first front of the associated wave, we 
estimated for each of CMEs 2-4 an average initial velocity in the range of 370 — 470 km s^^. 
We can speculate, that the kinematics shown in Fig. |2] (e) supports the idea of a strong 
deceleration of the faint, i.e. low-density, CMEs during their early evolution. 

In each of the four wave events, an increase in intensity in the central region of the 
source loop system (see Fig. I21^g)), constituting weak flares, is observed ~ 5 minutes before 
the first recorded wave front. Since the wave events occurred behind the solar limb, as 
viewed from Earth, no GOES X-ray data are available to obtain the flare class. As the 
determined wave kinematics allows us to back-extrapolate the start of the wave (see Fig. H] 
bottom panel), we can compare the timing of the flare commencement and the wave onsets 
(see also movie 2). Fig. H] (middle panel) shows the evolution of the total intensity of the 
central region in the AR. Four distinct intensity peaks are discernible, which coincide with 
the first observations of the four coronal waves. During these flaring phases we observe the 
disappearance of several loops in the 195A RR-images (see Fig. [2]^f)), whose northern branch 
is rooted close to the flaring part of the AR. The coronal waves are, however, launched from 
the opposite side and the shape of each first wave front exactly maps out the geometry of 
the southern loop system. 



4. Discussion and Conclusions 

Never before in the nearly 14 years of continuous research of large-scale coronal waves 
homologous EIT wave events were reported. We are the first to present a study of homol- 
ogous EIT waves, emerging from the same AR within a short period of time. They travel 
into the same direction and their fronts have similar shape and angular extent. They prop- 
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agate into a quiet Sun area, surrounded by ARs to the north and south and a large coronal 
hole close to the northern polar region (see movie 1). As is expected for nonlinear magne- 
tosonic waves, they do not penetrate into the se regions of increased Alfven velocity (see also 
Veronig et al.ll2008t iGopalswamy et al.ll2009[ ). In our study we compared for the first time 
different methods of deriving the wave kinematics, the (rather subjective) visual method and 
the (more objective) profile method. Both methods yield consistent results, i.e. the waves 
propagate at constant velocities ~ 220kms~^ for the weakest wave up to ~ 340kms~^ for 
the strongest event. Furthermore, we calculated the perturbation profiles to study the phys- 
ical characteristics and evolution of the disturbances. The strong initial steepening of the 
perturbation amplitudes confirms that these features are indeed shocks, albeit only weak 
shocks, since they peak at low intensity values A < 1.25. 

Assuming these coronal waves to constitute large-scale fast magnetosonic waves, the 
measured velocities lie well within the velocity range of 210-350 kms"^ for fast magnetosonic 
waves for quiet Sun conditions. In the MHD approach the quantities defining a shock wave 
in the solar corona are: Mms • • • shock magnetosonic Mach number, Xc . . . density jump at 
the shock front, ■(?... angle between shock front and magnetic field, /3c • • • plasm a-beta . The y 
are related by the Rankine-Hugoniot (RH) conditions for an oblique shock (c.f. |Priestlll982l ). 
Considering a perpendicular shock the RH relation reduces to 



'Xe(Xc + 5 + 5/3e 
(4-Xe)(2 + 7/3e 



with a polytropic index 7 of 5/3. Studies of IVrsnak et al.l (120021 ) indicate that (3 ~ 0.1 in 
the quiet Sun's low corona. Using the previously calculated density jumps Xc (cf. Sect. 13. 2[ 
Xc oc VA), we derive for the peak magnetosonic Mach numbers Mms= 1.06 (wave 1), 1.04 
(wave 2), 1.05 (wave 3), and 1.09 (wave 4). 

With these observations of homologous waves, we can for the first time perform a quan- 
titative analysis of the characteristic wave parameters without any limiting factors like 
changing or unknown quiet Sun background conditions. The top panel of Fig. [5] shows 
the calculated magnetosonic Mach number versus the propagation velocities Vc of the four 
waves, revealing a distinct correlation between the wave characteristics M^s and Vc with a 
correlation coefficient of ~ 0.99. Such correlation is expecte d for nonlin ear fast-mode 
magnetosonic waves because of the relation M^s = Vc/v^s (cf. iPriestI Il982l ). where Vc is 
the coronal wave speed and v^s = \A5"+~cf the fast magnetosonic speed [f a • • • Alfven 
velocity, Cg . . . sound speed] . The mean magnetosonic speed, derived by averaging the four 
ratios Vc/M^^, yields 250kms~^ and gives with = ISOkms"^ (T ~1.5 MK) a mean Alfven 
velocity of va = 175kms~^. Applying one- to five-fold Saito coronal density models for a 
propagation height of ~ 0.1-Rq, we derive for the low corona during solar minimum condi- 
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tions a magnetic field strength between 1.5 — 3.5 Gauss. 

The magnetosonic Mach numbers of all four waves peak at values less than 1.10 and 
drop off towards the linear regime (Mms ~ 1) as the waves expand. These small magne- 
tosonic Mach numbers together with the correlated wave velocities support the view that 
each observed coronal wave is a low-amplitude MHD fast-mode wave. Since M^s is even at 
the maximum close to the linear regime, this also implies that the fast-mode MHD waves 
are expected to experience only minor deceleration. Since Vc decreases proportional to Mmg, 
this implies a < 10% change in velocity for the waves under study. For the strongest wave 
4 this corresponds to a value of ~ 30 km s~^. As the error in velocity is of the same order 
of ~ 10% (~ ±30 km s~^ for wave 4), such weak deceleration is hidden in the measurement 
uncertainties. Finally, we note a clear correlation between the lags between successive waves, 
interpreted as the build-up times of the magnetic energy for the following wave event, and 
the intensity ///o, density N/Nq, magnetosonic Mach number Mms, and the propagation 
velocities Vc of waves 2 to 4 (cf. Fig. [5l bottom). 

In our study of the homologous wave events it was possible to compare the wave pulse 
characteristics and analyze correlations between the wave parameters, knowing that the 
four large-scale waves propagate in similar coronal conditions. This means that any found 
correlation is independent from the indeterminated but constant coronal conditions. Our 
results provide strong support that the observed large-scale coronal waves are indeed fast- 
mode magnetosonic waves, and additionally suggest a dependence of the wave parameters 
upon the build-up time of the magnetic energy. 

I.W.K., A.M. v., and N.M. acknowledge the Austrian Science Fund (FWF): P20867- 
N16. The European Community's Seventh Framework Programme (FP7/2007-2013) under 
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Fig. 1. — Median-filtered 10 minutes running ratio images recorded in EUVI-B 195 A showing 
the early evolution of four homologous waves on 2010 April 28-29 (axes in arcsec). Each 
column consists of two images recorded 20 minutes apart of one particular event. Upper 
panel: wave fronts at the time of the peak perturbation amplitudes [See accompanying 
movie 1]. Panel 7: yellow meridians define the sector, for which the wave kinematics was 
analyzed. 
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Fig. 2. — Associated features: (a)-(d) LASCO C2 running difference images of associated 
weak CMEs at the time of the maximum leading edge intensity [cf. movie 3]. (e) Kinematics 
of CMEs 2-4 in C2 (blue stars), plus estimated initial CME speed (green lines) by setting the 
moment of the first observed associated wave front as CME onset time. CME 1 was too faint 
to be measured, (f) Median-filtered EUVI-B 195 A running ratio image at 06:25 UT showing 
the early phase of wave 4 and the disappearance of an AR loop system (axes in arcsec). 
(g) Co-temporal direct image revealing an intensity increase of the inner loop system of AR 
11067 and the associated flare [cf. movie 2]. yellow square: 110" x 110" subfield used to 
analyze flare intensity variations. 
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Fig. 3. — Top: Kinematics of wave 4 determined by two methods plus error bars reflecting 
the diffusiveness of the wave fronts. Kinematics of the visually determined wave fronts in 
195 A represented by black circles, of the gaussian fitted perturbation profile by diamonds. A 
linear (black) and a quadratic (red) least square fit are overlaid. The grey area indicates the 
95% confidence interval of the linear fit. Bottom: evolution of the perturbation amplitude 
determined from EUVI-B 195 A intensity profiles in units of pre-event intensity Iq (compare 
Fig. m top panels). 



- 15 - 




22:00 00:00 02:00 04:00 06:00 08:00 

Time (UT) 



Fig. 4. — Top: 195 A peak intensity enhancements ('perturbation profiles') of tlie four waves 
togetlier witli Gaussian fits derived from RR images in a 60° sector (c.f. Fig. [H panel 7). 
Each vertical line indicates the distance of the visually tracked wave front, grey bars represent 
measured errors. Middle: variations of the flare intensity (in units of pre-event intensity; cf. 
Fig. Et^g), yellow square). Bottom: time distance diagram of all waves together with their 
linear fits. Green areas illustrate the 95% confidence bands for each linear fit. Note that the 
first wave fronts coincide with the flare peak intensity of the subfield. 
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Fig. 5. — Top: Correlation between magnetosonic Mach numbers Mms, derived from peak 
perturbation amplitudes, and linear propagation velocities Vc of the four coronal waves. 
Error bars in the Mach numbers are determined by the uncertainties in the gaussian fits to 
the perturbation amplitudes (see Fig. HI top panels). Bottom: Intensity amplitude (///o), 
density amplitude {N/Nq), Mach number (Mms), and wave velocity (fc) presented as a 
function of the time lags between successive waves. 



